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Background: A complete set is a minimum set of observables which allows one to determine the underlying reaction amplitudes unam- 
biguously. Pseudoscalar-meson photoproduction from the nucleon is characterized by four such amplitudes and complete sets involve 
single- and double-polarization observables. 

Purpose: Identify complete sets of observables, and study how measurements with finite error bars impact their potential to determine the 
reaction amplitudes unambiguously. 

Method: The authors provide arguments to employ the trans versity representation in order to determine the amplitudes in pseudoscalar-meson 
photoproduction. It is studied whether the amplitudes in the transversity basis for the yp —¥ K + A reaction can be estimated without 
ambiguity. To this end, data from the GRAAL collaboration and simulations from a realistic model are analyzed. 

Results: It is illustrated that the moduli of normalized transversity amplitudes can be determined from precise single-polarization data. Starting 
from simulations with achievable experimental resolution, it is quite likely to obtain imaginary solutions for the relative phases of the 
amplitudes. Also the real solutions face a discrete phase ambiguity which makes it impossible to obtain a statistically significant solution 
for the relative phases at realistic experimental conditions. 

Conclusions: Single polarization observables are effective in determining the moduli of the amplitudes in a transversity basis. Determining 
the relative phases of the amplitudes from double-polarization observables is far less evident. The availability of a complete set of 
observables does not allow one to unambiguously determine the reaction amplitudes with statistical significance. 

PACS numbers: ll.80.Cr, 13.60.Le, 24.10.-i, 25.20.Lj 



I. INTRODUCTION 

Pseudoscalar-meson photoproduction from the nucleon 
continues to be an invaluable source of information about the 
operation of quantum chromodynamics (QCD) in the non- 
perturbative regime [1-5]. Various classes of models, like 
constituent-quark approaches, have been developed to under- 
stand the structure and dynamics of hadrons in the low-energy 
regime of QCD. Experimental determination of the reaction 
amplitudes represents the most stringent test of those models 
and may open up a new chapter in the understanding of the 
energy eigenvalues and decay properties of hadron states. 

Quantum mechanics dictates that measurable quantities can 
be expressed as bilinear combinations of complex amplitudes. 
Pseudoscalar meson photoproduction involves only two kine- 
matical degrees of freedom, for example the energy of the 
incident photon and the scattering angle of the pseudoscalar 
meson. Since the (real) photon, the target, and the recoil- 
ing baryon have two spin degrees of freedom, eight possible 
amplitudes can be constructed. Due to angular momentum 
conservation, and depending on the adopted representation, 
half of these amplitudes either vanish identically or can be 
expressed in terms of the other four amplitudes. This leaves 
one with a set of four independent amplitudes. These ampli- 
tudes are complex functions (of the two kinematical variables) 
and therefore eight real functions are to be distinguished: four 
moduli and four phases. As quantum states are determined up 
to a constant phase factor, only the relative phases of the am- 
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plitudes can be extracted. This means that pseudoscalar pho- 
toproduction can be quantified in terms of seven real-valued 
functions of the two kinematical variables. Equivalently, at 
fixed kinematics seven real values suffice to determine all ob- 
servables. 

A set containing a minimum number of observables from 
which, at fixed kinematics, the seven real values can be de- 
termined unanimously is referred to as a complete set. In a 
seminal paper dating back to 1975, Barker, Donnachie, and 
Storrow argued that a complete set requires nine observables 
of a specific type [6]. In 1996, this was contested by Keaton 
and Workman [7] and by Wen-Tai Chiang and Tabakin [8]. 
The latter proved that eight well-chosen observables suffice to 
unambiguously determine the four moduli and three indepen- 
dent relative phases. 

In theory, a complete set of eight observables suffices to 
retrieve the generating seven variables. The need for an ad- 
ditional observable (a set of eight equations involving only 
seven variables) to uniquely determine the relative phases, is 
a reflection of the fact that the latter are linked to the observ- 
ables through nonlinear equations. A set of seven well-chosen 
observables generally yields multiple solutions for the phases. 
Therefore, an extra observable is required to mark the correct 
solution. In reality, however, observables have a finite preci- 
sion and deviate from the exact values. This compromises the 
solution of the phase ambiguity [9]. It then remains a question 
whether a set of eight observables is still sufficient to reach a 
situation of complete knowledge about the amplitudes. 

Thanks to recent technological advances in producing high- 
quality polarized beams and in developing polarized nu- 
cleon targets [10], it becomes possible to measure a suffi- 
ciently large amount of single- and double-polarization ob- 
servables in pion and kaon photoproduction. As a result, a 



2 



status of complete quantum mechanical information of pseu- 
doscalar meson photoproduction comes within reach. The 
self-analyzing character of the A is an enormous asset for 
achieving truly complete measurements for yp — > K + A and 
experimental efforts are under way [2]. For example, the 
CLAS collaboration at Jefferson Lab has many yp — > K + A 
polarization data in the pipeline. 

It should be stressed that in the quality and quantity of the 
experimental results, there is some kind of hierarchy. Indeed, 
double-polarization data is most often outnumbered by single- 
polarization results, which in their turn are outnumbered by 
the size of the differential cross section database. In view of 
this, the transversity representation of the amplitudes is a very 
promising one. Indeed, in this basis the single-polarization 
observables are linked to the squared moduli of the amplitudes 
by means of linear equations. Accordingly, the transversity 
basis occupies a central position in this work. 

The outline of the remainder of this paper is as follows. 
In Sec. II, the transversity amplitudes are introduced and all 
possible observables for pseudoscalar meson photoproduction 
are expressed in this basis. In Sec. Ill, Wen-Tai Chiang and 
Tabakin's formalism for solving complete sets is briefly re- 
viewed. This formalism is solid in the exact case. When ex- 
perimental uncertainty is involved, however, a consistency is- 
sue for a unique determination of the phases arises. A method 
to resolve this inconsistency is proposed. As a test of this 
method, in Sec. IV it is applied to simulations from a realis- 
tic model for the yp — > K + A reaction. The conclusions are 
summarized in Sec. V. In the Appendix, both the helicity and 
the CGLN expansions of the observables are considered, and 
a connection to the transversity basis is established. 



II. OBSERVABLES FOR PSEUDOSCALAR MESON 
PRODUCTION 

A. Reaction amplitudes 

The convention is adopted that the xz-plane coincides with 
the reaction plane and that the positive z-axis is along the 
direction of the photon's three-momentum. The two inde- 
pendent kinematic variables that will be considered through- 
out this work are the total energy of the reaction, or invari- 
ant mass, W and the scattering angle of the pseudoscalar 
meson in the center-of-mass frame c m .. As mentioned in 
Sec. I, pseudoscalar-meson photoproduction can be quanti- 
fied by four complex reaction amplitudes _M,(W,cos c . m .) 
(; = 1,2,3,4). The unpolarized differential cross section, for 
example, is given by 



d(j(W,cos0 cll 
dQ. 



= pY J \M i (W,co S e c . m .)\ 2 , (1) 

(=1 



where p is a kinematic factor. There are various equiva- 
lent representations for the M.\, all of which have a distinct 
kinematic factor p. The representation is determined by the 
choices made with regard to the quantization axis of the in- 
volved particles with a non-vanishing spin: the nucleon target, 
the incoming photon, and the recoiling baryon. 



The Dirac spinors for a particle with rest mass m, four- 
momentum p^ = (E.p), and a spin vector with polar angles 
(0,0) are defined as 



K±(^,6>,0) 



1 



yj2m{E - 



{E+m)I 2 
dp 



s±{0,<l>), (2) 



following the Bj0rken-Drell convention [11]. Here, 1% repre- 
sents the 2x2 identity matrix and a = (a x , a y , o z ) is the Pauli 
vector. Further, the Pauli spinors s + (0,0) ("spin up") and 
s_(0,0) ("spin down") are given by 



The following shorthand notations are introduced 



(3) 



|±), = «±(/Af,o), 
|±)_ v = «±(/Af,f) 
|±) z = «±(^,o,o). 



(4) 



Using Eq. (2) it can readily be verified that the |±) v and |±) 
can be expressed in terms of the |±) v as follows 



l±) 



(\-)y±\+)y), 



(5) 



The VW,(VK,cos c m .) are the amplitudes of the operator S^J^ 
for a fixed polarization of the initial and final state particles. 
The e£ is the photon's polarization four- vector and 7^ the 
transition current operator. The photon beam polarizations de- 
noted by "B = x" (+e x direction) and "B — y" (+e y direction) 
correspond with 



— 



(0,1,0,0), 
(0,0,1,0), 



and give rise to the following currents 



(6) 



Jx — &x ' 

Jy = EyJfl- 



(7) 



The photon beam polarizations denoted by "B — ±7r/4" 
(oblique polarization, £± K / 4 — (0, 1,±1,0)) and "B = ±" 

(circular polarization, = (0, 1, ±i, 0)) correspond with 



J±k/4 - ~^( Jx±J ^ ' 
l , 



(8) 
(9) 
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B. Observables in the transversity basis 



= ?(N 2 +N 2 



(14) 



The so-called transversity amplitudes bi express the Mj in an( j 
terms of the spinors |±) y (quantization axis perpendicular to 
the reaction plane) and the linear photon polarizations J x and 
J y , i.e. 

b\ =v( + |A-|+>y, 
b2=y(-\Jy\-)y, 
b 3 = y (+\Jx\-)y, 

b 4 =y(-\J x \+)y. (10) 

In these definitions, the bra (ket) refers to the recoil (target). 
The differential cross section for a given beam £>, target T, 
and recoil polarization 7Z is denoted as 



.(b,t,u) 



di2 



An unpolarized state is denoted by "0". For example, T = 
denotes an unpolarized target and in computing the cross sec- 
tion (11) for T = 0, an averaging over both target polarizations 
is implicitly assumed. 

An asymmetry A can generally be expressed as 



A = 



g{B\,T\,K X ) _ g{B 2 ,T 2 ,TL 2 ) 

g(Bi,ri,TC 1 ) + g(B2,r 2 ,7j 2 )- 



(12) 



A single asymmetry comprises one polarized and two unpolar- 
ized states. Hence, there are three possible single-polarization 
observables, namely the beam asymmetry E (B ^ 0), the target 
asymmetry T (T ^ 0), and the recoil asymmetry P (1Z ^ 0). 
The explicit definitions of these three single asymmetries can 
be found in Table I. A double asymmetry, involves two polar- 
ized and one unpolarized state. There are three types of double 
asymmetries: the target-recoil asymmetries (By = B 2 — 0), the 
beam-recoil asymmetries (71 = Ti = 0), and the beam-target 
asymmetries (IZi = 1Z 2 — 0). The definitions for the various 
double asymmetries are contained in Table I. 

The aim of the current section is to represent the single and 
double asymmetries in the transversity basis. As representa- 
tive examples, the transversity expansions of the single asym- 
metry E and of the double asymmetry C x are derived. The 
beam asymmetry E is defined as 



gCvAO) _ g(*,0,0) 

g(y,o,o) + g(*,o,or 



(13) 



Using 



.(y,0,0) 



)4-gb>>+y-y) 



.q(y,-y,+y) _|_ gty>-y>-y. 



y(+\Jy\+)y 

v v ' 



2 +\y(-\Jr\+)y\ 2 



=0 



+ \y{+\Jy\-)y\ 2 +\y{-\Jy\-)y\ 2 ), 



(x : 0,0) 



+y) 4- A x >+y -y) 



. g( x -y>+y) _|_ g( x -y-y) 











{\y(+\J 


x\+)y\ 2 ^ 


-\A-\J* 


\+)yf 











=0 



+ \y{+\Jx\-)y\ 2 +\y{-\Jx\-)y\ 2 ), 



2 



^(\b3\ 2 + \b 4 \ 



(11) one obtains 



E = 



Ib^ + M^ib^-M 2 



I^l| 2 + |fe2| 2 + |^| 2 + |fe4| 2 ' 

At this point the normalized transversity amplitudes a,- 

bi 

°' V\bi\ 2 + \b2\ 2 + \b3\ 2 + \b4\ 21 



(15) 



(16) 



(17) 



are introduced. Upon using the notation r,- = the normal- 
ization condition of the a, reads 



r\ + rl+rl + rl = l, 



(18) 



which means that there are only three independent moduli, or 
six real values to be determined. For the bi, this is respectively 
four and seven. The \bi\ can be obtained from the r,, combined 
with the magnitude of the unpolarized differential cross sec- 
tion. Using the r ; -, the beam asymmetry E can be expressed 
as 



v 2,2 2 2 

?- = r l +r 2 -r i -r A . 



(19) 



Similar calculations as those for E yield the expressions for T 
and P, contained in Table I. 

Next the beam-recoil asymmetry C x 



C, 



g(+,0,+x) _ g(+,0,-x) 
g( + fi,+x) + g( + .0.-x) > 



(20) 



is considered. The first term in the numerator can be expressed 
as 



-.0, 



"4( 



g(+,+y,+X) + g( 



~-b 2 



Using Eq. (9), one obtains 

g {+ '°' +X) = ^(\x(+\Jx\+)y 
+ Ix<+|A|— )vf 



\Jy\ 



(21) 



\Jy\ 



(22) 
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TABLE I: The expressions for the single and double asymmetries in the normalized transversity basis. Both the expressions obtained in this 
work and those found in the literature [13] are listed. The convention for the beam-target asymmetry E is adopted from Ref. [12] instead of 
Ref. [13]. 





{B h Ti,Ki) 


(B2,T 2 ,n 2 ) 


Transversity representation 

This work Literature [13] 


£ 


(y,o,o) 


(x,0,0) 






i- 2 — r 2 
r 3 r 4 


T 


(0,+y,0) 


(0,-y,0) 






r 2 + r 2 


P 


(0,0, +y) 


(0,0, -y) 






r\-r\ 


c x 


i+A+x) 


(4-,0,-x) 


— 2Im(aifl4 


-1-0203) 


— 21m (a 1 — 02^3 ) 


c z 


(+,0,+z) 


(4-,0,-z) 


+2Re(aiaJ 


-0203) 


+2Re (a 1 + 0203 ) 


o, 


(+f,o,+*) 


(+f,0,-x) 


i on„/„ _* 

-\-lKe{a\a^ 




+2Re(aici4 — 0203) 


o z 


(+f,o,+z) 


(+f,0, -z) 


+2Im(aifl4 


-a 2 fl3) 


+2Im (a 1 a| + fl2 a 3 ) 




f i ^ n^ 




+2Re(aia^ 


- 0204) 


— 2Re(aiflj +02^4) 


F 


(+,+*> o) 


(+.-*, o) 


— 2Im(aifl3 


f a 2 a\) 


+2Im(aia? — 02^4) 


G 


(+f,+z,0) 


(+f -z,0) 


— 2Im(aia^ 


- a 2 a\) 


+2Im(aifl3 40204) 


H 


(+!,+*, o) 


(4-f,-*,0) 


+2Re(a ia ; 


+■ 02^4) 


— 2Re(ai«3 —0204) 


T, 


(0, +*,+*) 


(0, +x, —x) 


+2Re(ai02 


fa 3 a!) 


+2Re(«ia2 — 03O4) 


T z 


(0,+x,+z) 


(0,+x,-z) 


+2Im(aifl2 


-1-0304) 


+2Im(aia2 — 03O4) 


L x 


(0,4-2,4*) 


(0,+z,-x) 


— 2Im(ai<a!2 


- 0304) 


— 2Im(aifl2 + 03O4) 




(0,+z,+z) 


(0,+z,-z) 


+2Re(aia| 


-0304) 


+2Re(ai«2 + a 3o|) 



As it is unfrequently used in the field, the transversity rep- 
resentation for the observables is not contained in the exten- 
sive recent review by Sandorfi et al. [12]. However, when it 
comes to extracting the invariant amplitudes from complete 
measurements, it will turn out that the transversity represen- 
tation is highly beneficial. This will the subject of Sec. IV B. 
The transversity expressions for the double asymmetries ob- 
tained in this work, are not consistent with those listed in lit- 
erature [6, 9, 13]. 



By inspecting Table I it is clear that the substitution 03 — > 
—03 makes the expressions derived in this work consistent 
with those contained in Ref. [13]. To the knowledge of the au- 
thors, the expressions for the asymmetries in the transversity 
basis were first published by Barker, Donnachie, and Storrow 
in 1975 [6]. Slightly different expressions are contained in a 
1990 paper by Adelseck and Saghai [13]. However, explicit 
derivations are not contained in either of both publications. 
In Appendices 1 and 2, the asymmetries are expanded in the 
helicity and the CGLN basis, respectively. The helicity rep- 
resentation is presented in Table VIII and is consistent with 
the one obtained by Fasano et al. [14]. On the other hand, 
the CGLN expansion was found to coincide with the one cal- 
culated by Sandorfi et al. [12]. Interestingly enough, both 
the helicity and CGLN representations were derived from the 
transversity expressions obtained in this work. Starting from 
the transversity expressions listed in the last column of Table I 
the helicity and CGLN expansions available in the literature 
are not retrieved. 



Then, Eq. (5) is employed to transform the x (+\ spinors into 
v (±| spinors, i.e. 



(\y(+\Jx\+)y 


+y(-\Jx\+)y 


= 


=b A 


iy(+\Jy\+)y- 


Hy(-\Jy\+)y 


= bl 


= 




-yHJxhh 


=63 


= 


iy(+\Jy\-)y- 


S v ' 


=0 


=b 2 


(|^l+£4| 2 -+ 


-\ib 2 + b 3 \ 2 ). 



(23) 

An analogous calculation for g(+>° ~*) yields 

g (+,o,- x ) = I (\ ibl -b 4 \ 2 + \ib 2 ~b 3 \ 2 ) . (24) 

Inserting expressions (23) and (24) in the definition (20) leads 
to 

i{bib* 4 -b* l b4 + b 2 b* 3 -blb 3 ) 
x ~ N 2 4-N 2 4-|fc 3 | 2 + N 2 ' 
= — 2Im(aia4 + 0203). (25) 

Analogous derivations for the remaining beam-recoil 
(C z ,O x ,O x ), the beam-target {E,F, G,H), and the target-recoil 
asymmetries (T X ,T Z ,L X ,L Z ) yield the expressions listed in 
Table I. 
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III. INFERRING THE TRANSVERSITY AMPLITUDES 
FROM THE ASYMMETRIES 

A. The moduli 

In Sec. IIB the notation r, = |a,| was introduced. Hence, 
the (normalized) transversity amplitudes can be expressed as 

a, }■;<■". (26) 

Here, a, represents the phase of the transversity amplitude a,-. 
From the expressions for E, T and P in Table I and from the 
normalization condition (18), one can readily solve for the 
moduli of the amplitudes in terms of the single asymmetries, 
to obtain 

1 , 

~ P, 



r\ = 



>~2 



'"3 



>~4 



T.-T-P, 



= -Vl-H-T+P, 



Vl-Z + T-P. 



(27) 



Hence, a measurement of E, T and P at fixed kinematics 
(W, cos0 c m .) yields the moduli of the transversity amplitudes 

n- 

In Sec. IIB, it was mentioned that only three of the four 
moduli are independent. Any of the four combinations of 
three independent moduli can be chosen without violating 
generality. As will be explained in Sec. IV C 2, however, it 
is not beneficial to eliminate one of the moduli. 



B. The relative phases 

1. Independent and nonindependent relative phases 

From the double asymmetries one can extract the relative 
phases a, ; - = a, — a,, given that the moduli r, are known. 
There are six possible combinations for the relative phases, 
namely an, «13, OCu, Q-23, &24, and 0:34, all of which can be 
extracted from the double asymmetries (see Table I). How- 
ever, only three of these are independent variables; the re- 
maining three can be expressed as linear combinations of the 
three independent relative phases. For example, with OC4 as the 
reference phase, (ai4, 0:24, 0:34) are the independent phases 
and (0!i2, 0!i3, 0:23) are the nonindependent ones 

«12 = 0:14-0:24, 

^3 = 0!l4-0!34, 

«23 = «24-«34- (28) 

Alternatively, one can express the phases relative to (X\, OC2, or 
a.3 . From now on, 0:4 will serve as the reference phase and the 
independent phases are denoted as <5, = 05; — 0:4 (i = 1,2,3). 
The nonindependent phases are labeled as A,y = 8, — 8j (i ^ 
j). Table II shows the correspondence between the indepen- 
dent and the nonindependent phases for the four distinct ref- 
erence phases. 



2. Complete sets 

In Ref. [8], Wen-Tai Chiang and Tabakin proved that a spe- 
cific set of four double asymmetries is sufficient to extract two 
of the independent phases, 5, and 8j, and two of the nonin- 
dependent phases, A# and A,^ (z ^ ^ k). The remaining 
independent phase 5^ can be constructed as 8k = 5, — A,* or 
8^ = 8j — A jit. Hence, a complete set, consisting of three sin- 
gle and four double asymmetries, determines the four moduli 
and the three independent relative phases of the transversity 
amplitudes a,-. 

The collection of complete sets, which is listed in tables 
III— VIII of Ref. [8], can be divided into two categories. An 
example of a complete set of the first kind is {C x ,O x ,E ,F} 
(the single asymmetries E, T, and P are implicitly assumed to 
be included as they are a vital part of any complete set in the 
transversity basis). From Table I one obtains 



and 



ri^4 sin 81 + nr?, sinA23 : 
r\ 7-4 cos 8\ + r2f"3 cos A23 = 

ri^cosAi^ — r2f"4Cos 82 
rir3sinAi3 + ^^1x182 



Cx 

2 ' 

Ox 

2 ' 



(29) 



E 
2' 



F 
2' 



(30) 



Since Eq. (29) contains both a sine and a cosine of the un- 
knowns, there are generally two solutions for {5i,A23}. The 
same reasoning applies to E and F, and hence there are four 
possible solutions for {81 , 82, A13, A23}. Yet, only one of these 
solutions will satisfy the (trivial) relation 



<5i+A23-<5 2 -Ai3=0, 



(31) 



and that specific solution is the actual solution. 

If, for example, one replaces F with H in Eq. (30), a com- 
plete set of the second kind is obtained. Table I yields 



nr3COsAi3 +r2?"4COS(52 



H 



(32) 



Since {E,H} contain only the cosine of the unknowns, there 
are four possible solutions for {(>2,Ai3}. Hence, an eight- 
fold ambiguity for {5i,(>2,Ai3,A23} emerges. Only one of 
the solutions, however, obeys the constraint (31). Similarly, 
the complete set {C x ,O x ,F, G} is subject to an eightfold phase 
ambiguity. 



C. Introducing experimental error 

As indicated in Sec. IIIB2, a complete set provides ac- 
cess to two independent and two nonindependent phases, i.e. 
the set {8i, 8j, A,*, A^} (z ^ j ^ k). There are two ways to 
calculate the third independent phase 8k'- 8k — 5, — A,* or 
8k = 8j — Ajk- For infinite experimental resolution both ex- 
pressions for 8k are equivalent. When experimental error is 
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TABLE II: Correspondence between the independent (5,) and nonin- 
dependent (Ay) phases for the reference phase <Xt\ and the indepen- 
dent (8 l at ) and nonindependent (Ay = 8 i at — 5 * = ctj — CCj) phases 
for the reference phase (k= 1,2,3). 

«4 ^2 #3 A12 Alj A 23 

a, -8 4 a ' A 24 A 34 -5 2 a ' -5 3 a ' A 23 
«2 A 14 -V 2 A 34 «/* A B -5 3 " 2 
a 3 A 14 A 24 -5 4 " 3 A12 V 3 g 2 " 3 



introduced, the estimated values {5,-, 8j, Ai k) A»} deviate from 
their actual values 5y, A,*, A^} (estimators are marked 

with a " " ") and the two expressions for 8^ yield a differ- 
ent result in general. Since there is no prior preference, an 
equal weight can be assigned to both estimators 

8 k = i(8 i -A ik ) + ^(8j-A jk ), (33) 

thereby resolving the 8^ ambiguity. However, this is not 
the only problem that emerges from introducing experimen- 
tal error. Consider again the complete set {C x ,Ojc,E,F} of 
Eqs. (29)-(30). For this set of observables the 5, read 

% 

< 82, (34) 
^ = i(5 1 -Ai 3 ) + i(^-A 23 ). 

Suppose that, for example, a,\ was assigned as the reference 
phase, instead of OC4. Then, {C x ,O x ,E,F} provides access to 

{S," 1 , 5 4 a ' , A23, A24}, as can be inferred from Table II. The 
corresponding independent set of phases reads 

p 2 "> = 1(^+^3) + ^+^), 

<S 3 a> , (35) 

From (35) one can estimate the 5, through {8{ = — 8 4 a ' , 8 2 ' = 
8 2 " 1 - 8 4 ai , 5/ = 8" 1 - 8 4 l } and Table II. This yields 

{8{ = 8 U 

i^=i%+i(5 1+ A 23 -A 13 ), (36) 
[s 3 ' = Si-A 13 . 

The set (36) is not consistent with the set (34). Indeed, the 
estimates for the independent phases (there are four possible 
sets) depends on the choice of reference phase. 

One would like to have a consistent set of estimators for 
the independent phases 5. ' (i 7^ j) and S" 4 = 5, (i — 1,2,3). 
The notation 8^' is adopted for the consistent estimators. The 



aforementioned reference-phase ambiguity can be resolved by 
imposing that 

8" J = 8 i ak -8 J CCk . (37) 

Since this requirement is linear in the independent phases, 
the consistent estimators are also linear in {5 ( a ' , 5^*' , A#, A^} 
(z 7^ j 7^ ky^ I). The most general expression for a consistent 
set of estimated phases, reads 

f S> = ei5,°» + (1 - ci) (8? + A ik - A Jk ) , 

I Sp =c 1 8 j a ' + (l-c 1 )(8 i a ' + Ajk - A,,) , (38) 

(8«=c 2 (8«-A ik ) + (l-c 2 )(S?-Aj k ). 

Here, three constraints have been imposed. Firstly, the 8"'. ^ 
should yield the exact set 8"'. ^ for vanishing experimental 

error bars. Secondly, the expressions for 5 ( a ' and 8"' should 
be equal for i O j. This leaves one with two unknowns c\ and 
c 2 . Thirdly, c\ and c 2 cannot depend on a/, which is implied 
by Eq. (37). By applying Eq. (37) on expressions (38) and by 
using Table II, one readily finds that ci = | and c 2 = \. There- 
fore, the consistent set of estimators for the phases, reads 

(A a ' = lSr' + l(8j a '+A ik -Aj k ), 

hp = l8 J ai + l(8 i ai +A jk -A ik ), (39) 

[^ = I(^-A ifc ) + I(5/'-A ; ,). 

The estimates 5,-'s (i = 1,2,3) for the independent phases are 
now insensitive to the choices made with regard to the refer- 
ence phase. 

IV. RESULTS 

A. The angular and energy dependence of the transversity 
amplitudes 

A complete measurement comprises a minimal set of asym- 
metries from which the accessible parameters of the normal- 
ized transversity amplitudes a, can be estimated. These pa- 
rameters include three independent moduli r, and three inde- 
pendent relative phases 5,-. In Sec. Ill, it was shown how the 
moduli can be obtained from the single asymmetries and how 
a set of three independent phases can be estimated in a consis- 
tent way from a complete measurement which also involves 
double asymmetries. 

Figures 1 and 2 show the (W, cos c .m.) dependence of 
the r,- and the 5, for the yp — > K + A reaction as predicted 
by a realistic model, namely the latest version of the Regge- 
plus-Resonance (RPR) model, i.e. RPR-2011 [15-17]. This 
model has a Reggeized f-channel background and includes a 
total of 8 i-channel resonances, namely Si i( 1535), Si i( 1650), 
Fi 5 (1680), A 3 (1720), At (1900), P 13 (1900), £>i 3 (1900), and 
^15(2000). The most apparent feature of the RPR-2011 pre- 
dictions for the moduli and phases are the strong variations 
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FIG. 1 : (Color online) The energy and angular dependence of the moduli r,- of the normalized transversity amplitudes for the yp — > K + A 
reaction. The calculations are performed with the RPR-201 1 model (AW & 2.38 MeV, Acos c . m . & 8.33 x 10~ 3 ). 



with energy W at backward scattering angles, due to s-channel 
resonances. The smooth energy dependence of the moduli 
and phases at very forward scattering angles reflects the im- 
portant f -channel background contributions of yp — > K + A. At 
forward kaon angles, where most of the strength resides, the 
RPR-201 1 model predicts a dominant role for r%. 



B. The merits of the transversity basis 

In Sec. Ill, it was shown how the transversity amplitudes 
can be inferred from a complete set of measured observables. 
The first and essential step comprised the determination of 
the moduli from the three single asymmetries. Indeed, as 
became clear, the extraction of the 5,- requires prior knowl- 
edge about the moduli r,-. Luckily enough, single asymmetries 
are more easily obtained experimentally than double asym- 
metries. Consequently, the single asymmetries database gen- 
erally have better statistics compared to double polarization 
observables. The published database for e.g. yp — > K + A in- 
cludes 2260 single asymmetries (£: 178, T: 69, P: 2013), 
in contrast to only 456 double asymmetries (C x : 162, C z : 
162, O x r. 66, O z r. 66) [17]. The beam-recoil asymmetries 
{O x i,O z i} are related to {O x ,O z } through Eq. (A.9). In the 
transversity representation the parameters that are most easily 
extracted from the data, namely the moduli, are related to that 



class of asymmetries that are more readily measured, namely 
the single asymmetries. This is not the case for the helicity 
basis, where the moduli are related to the set of double asym- 
metries {C z i ,E ,L z i}. It is important to stress that determin- 
ing invariant amplitudes at some kinematical point requires 
knowledge of a complete set at the same kinematics. To date, 
however, there is not a single pseudoscalar photoproduction 
reaction for which a complete dataset has been published. 

For the yp — > K + A reaction, the GRAAL collaboration has 
measured {£, T,P} at 66 kinematical points [18, 19]. These 
data cover 1.65 < W < 1.91 GeV (bins of AW w 50 MeV) and 
—0.81 < cos 9 c .m. ;$ 0.86 (bins of Acos c m . w 0.3), and can 
be used to estimate the moduli r, through Eqs. (27). Figure 
3 shows these extracted moduli at three cos 9 C m . intervals (33 
of the 66 kinematical points for which data is available) along 
with the corresponding predictions from RPR-201 1. It is ob- 
served that at a few kinematical points, some of the r, cannot 
be extracted from the data. This occurs whenever the mea- 
sured set {£, T,P} deviates too strongly from the 'exact' set 
so that one or more arguments of the square roots (27) become 
negative. For the GRAAL data, only 17 of the 66 x 4 = 264 
f; values are imaginary (r\\ 12, r^- 3, r^. 2). Overall, the 
RPR-201 1 model offers a fair description of the energy and 
angular dependence of the extracted r,- except for the Pi and 7\ 
at 0.808 < cos c m . < 0.861. 

Obviously, the GRAAL bin width of AW k, 50 MeV suf- 
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FIG. 2: (Color online) The energy and angular dependence of the 
independent phases 5, of the normalized transversity amplitudes of 
the yp — > K + A reaction. The calculations are performed with the 
RPR-2011 model (AW w 1.43 MeV, Acos6 c . m . = 5 x 1(T 3 ). 



{L,T,P, O x /,O z r} cannot be used to extract the moduli of the 
helicity amplitudes. The latter would require experimental 
data for {C z i,E,L z /} at the same kinematics, which is not 
available to this day. Clearly, given the current status of the 
experimental yp — > K + A program, the transversity represen- 
tation offers the best perspectives to learn about the reaction 
amplitudes. 

The extracted information about the energy and angular de- 
pendence of the amplitudes r,(W, cos cm .) of Fig. 3 is com- 
plementary to what could be obtained about the partial waves 
in a so-called truncated partial wave analysis (TPWA) [21- 
23]. In TPWA one aims at extracting information about the 
energy dependence of the partial waves from the experimen- 
tal data. The partial waves serve as expansion parameters for 
the angular dependence of the observables. The amount of 
partial waves which can be included (usually denoted by Z max ) 
in the fits, depends on the quality and quantity of the measured 
cos C m . dependence of the observables. It is anticipated that 
in a reaction channel like yp — > K + A with substantial back- 
ground contributions, / max cannot be truncated to small values. 
In the TPWA approach the underlying dynamics is parameter- 
ized in terms of functions which depend on W, whereas in an 
amplitude analysis as presented here, one attempts to map the 
full (W, cos Oc.m.) dependence. In the forthcoming section, it 
will be shown that the availability of complete sets of eight ob- 
servables with realistic error bars at given (W, cos c .m.) does 
not guarantee that one can retrieve the underlying transver- 
sity amplitudes. In a recent analysis of yN — > nN simulated 
data [21] it was shown that the availability of the cos 6 c m , de- 
pendence of six observables at given W yields unique partial- 
wave solutions. 



C. Extracting the transversity amplitudes from RPR-2011 
simulations 

In this section the completeness of complete sets is investi- 
gated. The analysis is performed for the representative com- 
plete set {C x ,O x ,E ,F} of the first kind. The goal is to deter- 
mine the moduli % and the consistent independent phases 8j of 
the transversity amplitudes from simulated observables with 
finite experimental resolution generated with the RPR-2011 
model. It is investigated to what extent the retrieved ampli- 
tudes comply with the input amplitudes from the simulations. 



fices to map the energy dependence of the r,. The most ex- 
tensive yp —> K + A data set to date is due to the CLAS col- 
laboration at Jefferson Lab [20]. The data cover 1.62 < W < 
2.84 GeV (bins of AW = 10 MeV) and -0.85 < cos c m . < 
0.95 (bins of Acos 9 C m . = 0.1) and include data for ^ and the 
recoil polarization P. 

In addition to {£, T,P}, GRAAL provides data for 
{O x i,O z i} at the same kinematical points. Since a complete 
set of asymmetries is required to extract the 5,, these beam- 
recoil data cannot be employed in order to gain additional in- 
formation about the transversity amplitudes. 

As is seen from Table VIII, the GRAAL data for 



1. Strategy 

A measured asymmetry is simulated by generating a fixed 
number of events from a Gaussian distribution. The mean of 
this distribution is the RPR-201 1 prediction for a certain value 
of W and cos c m . The standard deviation is a specific exper- 
imental resolution (7 e xp- The mean and standard deviation of 
the generated events determine the value and the error of the 
simulated data point. If the error of the simulated data point is 
smaller than the specified ff e xp, the former is rejected. 

From a simulated set {C x ,O x ,E ,F}, one can obtain the r; 

and the 5, by means of Eqs. (27) and (39). In order to esti- 
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FIG. 3: (Color online) The energy dependence of the moduli r; of the normalized transversity amplitudes for the yp — > K + A reaction. The 
data are extracted from the GRAAL results for the single-polarization observables reported in Refs. [18, 19]. The dots are the bin-centered 
RPR-201 1 predictions. 



mate the error on the r, and the 5,, standard error propagation 
is applied. Since seven asymmetries are required to estimate 
the six parameters of the normalized transversity amplitudes, 
the correlation between the asymmetries has to be taken into 
account upon estimating errors. The squared error of a certain 
function / of the asymmetries is calculated as 



- 2 (/)=zflfV- 2 ^) 



i,j • j 

¥J 

with A u Aj G {L,T,P,C x ,O x ,E,F} and a(Ai,Aj) the covari- 
ance between A, and A;. 



2. Analysis 

As indicated in Sec. IV B, the extraction of the r, from a 
measured set of single asymmetries through Eqs. (27) is rather 
straightforward. In some isolated situations at least one of the 
retrieved r, became imaginary. The moduli estimates, how- 
ever, are required to extract the relevant phases from a com- 
plete set, as is seen from Eqs. (29), (30), and (32). Conse- 
quently, an imaginary estimate for one of the moduli results 



TABLE III: The model's values and two simulated datasets for the 
complete set {C x ,O x ,E,F} at W = 1700 MeV and cos0 c . m . = -0.5 



with a. 



exp 



0.05. 





Model 


Simulation A 


Simulation B 


z 


0.0489 


0.0543±0.0519 


0.0458 ±0.0533 


T 


-0.6843 


-0.6905 ±0.05 11 


-0.6685 ±0.0524 


P 


0.0056 


0.0111 ±0.0552 


0.0892 ±0.0555 


c x 


-0.4808 


-0.4809 ±0.05 15 


-0.4728 ±0.0588 


o x 


-0.5989 


-0.5971 ±0.0519 


-0.6149±0.0567 


E 


0.5672 


0.5655 ±0.0505 


0.5695 ±0.0564 


F 


-0.4525 


-0.4411 ±0.0519 


-0.4486 ±0.0503 



in imaginary estimates for the phases, no matter the achieved 
precision of the double-asymmetry observables. Through the 
normalization condition (18) only three of the four moduli are 
independent. One imaginary modulus estimate, be it an inde- 
pendent or a nonindependent one, is sufficient to jeopardize 
the phase analysis. For this reason, the four moduli are treated 
on equal grounds and a selection of a set of three independent 
moduli has been avoided all along. 

Estimating the 5,- is far less straightforward than estimat- 
ing the r,-. Table III lists two simulated measurements of the 
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complete set {C x ,O x ,E,F} at W = 1700 MeV and cos c . m . = 
—0.5, performed at an input resolution of (7 e xp = 0.05. The 
mean value and standard deviation of the asymmetries were 
calculated from a sample of 50 events. Upon inspecting Ta- 
ble III, both simulations are seen to be qualitatively equiv- 
alent. This is also reflected by the accurate moduli esti- 
mates from both simulations, which are listed in Table IV. 
In Sec. IIIB 2 it was mentioned that four distinct solutions for 
{5i,(>2,Ai3,A23} exist for the concerned complete set. Only 
one satisfies the constraint (31) in the exact case, though. 
When experimental error is involved, however, none of the 
four solutions satisfies the constraint. In Table V, the outcome 
of constraint (3 1) is listed for each of the four solutions to both 
datasets. As the constraint is a sum of phases and its eval- 
uation in the exact case yields a zero value, the mean values 
listed in Table V were rotated to ] — n, +n] . The table also lists 
by how many standard deviations (denoted by n a ) the mean 
deviates from zero. The solution with the lowest n a has the 
highest likelihood. In principle a solution can be excluded if 
it is statistically insignificant, commonly quantified by a small 

confidence interval CI = 1 — erf (n a / \flj . From the CI val- 
ues of dataset B it follows that solution 3 and solution 4 can be 
excluded with significance. For dataset A, these solutions are 
also the least significant ones. Of the two remaining solutions, 
dataset A predicts that solution 1 is the real solution (99.4% 
Cl). According to dataset B, however, solution 2 is the most 
likely one, though with a smaller confidence interval (79.9% 
Cl). This example illustrates that ambiguities remain for (the- 
oretically) complete sets when experimental uncertainties are 
taken into account. 

In Table VI, the model's phases are compared with the cor- 
responding estimates from solution 1 of dataset A and solu- 
tion 2 of dataset B. From this table, it is seen that solution 1 is 
the actual solution. This finding could have been anticipated 
by comparing CI values (99.4% versus 79.9%). However, a 
real experiment only yields a single dataset at a certain kine- 
matical point. If an experiment were to yield dataset A, one 
could not exclude solutions 3 and 4 with significance and one 
could only conclude that solution 1 is most likely the real so- 
lution with 99.4% CI. If dataset B were the result of a real 
experiment, then solutions 3 and 4 could be excluded with 
significance and one could state that solution 2 is most likely 
the actual solution with a CI value of only 79.9%, though. The 
confidence interval on the constraint (3 1 ) is the only parameter 
able to distinguish between the four solutions. In some cases, 



TABLE V: The constraint of Eq. (31) for the four {Si , &i, A 13 ,A 23 } 
solutions extracted from the simulated data listed in Table III. The 
number of standard deviations, by which the constraint differs from 
zero, is denoted by n a - The CI are the confidence intervals corre- 
sponding with the n a . 



Solution 


Simulation A 

Constraint n a 


CI (%) 


1 


0.029 ±3.866 


0.007 


99.4 


2 


-0.200±3.316 


0.060 


95.2 


3 


-2.513±3.316 


0.758 


44.9 


4 


-2.741 ±3.865 


0.709 


47.8 



Solution 


Simulation B 

Constraint n a 


CI (%) 


1 


0.592± 1.261 


0.467 


63.8 


2 


-0.229 ±0.898 


0.255 


79.9 


3 


-2.393 ±0.903 


2.651 


0.8 


4 


3.069 ±1.264 


2.428 


1.5 



however, the solution with the highest CI does not deliver the 
correct solution for the phases. 

In the beginning of this section, it was argued that an imag- 
inary estimate for the r, leads to imaginary 5;. Even with real 
r,-, however, it cannot be excluded that imaginary 5, are re- 
trieved. Indeed, from Eqs. (29) and (30) the phases are de- 
termined by solving a quadratic equation in the sine or cosine 
of the phase. For certain combinations of the measured dou- 
ble asymetries the discriminant of this quadratic equation be- 
comes negative. In such case, the particular dataset can solely 
be used to estimate the moduli. In what follows the authors 
wish to quantify the insolvability of a complete set, i.e. how 
frequently it occurs that imaginary or incorrect solutions are 
obtained for the trans versity amplitudes. 

Figure 4 shows the insolvability of the transversity am- 
plitudes as a function of the input experimental resolution 
(7 exp . Simulations for the complete sets {C x ,O x ,E,F} and 
{C x , O x , G, H} are analyzed at four different kinematics. Each 
point in Fig. 4 results from an ensemble of 200 simulated 
datasets constructed from samples of 50 events for each of 
the asymmetries. The insolvability rj is the fraction of un- 



TABLE IV: The model's values for the % at W = 1700 MeV and 
cos f c .m. = — 0.5, and the corresponding estimates from dataset A 
and dataset B. 





Model 


Simulation A 


Simulation B 


n 


0.304 


0.306 ±0.037 


0.311 ±0.037 


h 


0.657 


0.658±0.017 


0.653 ±0.018 


h 


0.640 


0.642±0.018 


0.639±0.018 




0.256 


0.248 ±0.046 


0.263 ±0.044 



TABLE VI: The model's values for the 5, at W = 1700 MeV and 
cos0 c m . = —0.5, along with the estimates resulting from the most 
significant solution of dataset A (solution 1) and dataset B (solution 
2). 







Model 


Simulation A 


Simulation B 


Si 


0.487 


0.424 ±0.801 


0.704 ±0.478 


Si 


2.458 


2.424 ±0.992 


2.637 ±0.348 


5 3 


6.106 


6.090 ±0.235 


3.046 ±0.234 
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FIG. 4: (Color online) The insolvability of the transversity amplitudes at four different kinematical points. The circles and triangles are 
simulated data generated from the complete sets {C x ,O x ,E,F} and {C x ,O x ,G,H}. The filled and unfilled symbols correspond with 77 = 

^imaginary ^incorrect <*nd 77[ ma g[ nal -y ■ 



successful simulated datasets. The tj is defined as the sum of 
the fraction of (two or four) imaginary solutions (^imaginary) 
and the fraction of incorrect solutions (i7incorrect)- An incorrect 
solution is a solution with the highest CI value that does not 
correspond with the model's value in the limit ff eX p — > 0. The 
insolvability 7} is binomially distributed and hence the corre- 
sponding error is calculated as 

where N represents the number of simulated datasets. It is 
seen that the insolvability decreases with decreasing a exp , as 
is expected. Further, the I7i ma gmary has the largest contribution 
and i7i n conect vanishes more rapidly with ff exp than I7i ma ginary- 
In some situations, the insolvability only vanishes for very 
challenging experimental resolutions, and it is seen that this 
behavior depends on the complete set in question. 

Figure 5 shows the {C x ,O x ,E,F} insolvabilities 7} and 
^incorrect as a function of W and cos c m for dexp = 0.1 and 
cj exp = 0.01. For each of the 375 x 241 = 90375 kinematical 
points, 1000 simulated datasets were generated (from sam- 
ples of 50 events for all asymmetries) from which the av- 
erage T7i m aginary an d ^incorrect are calculated. It is seen that 
for ffexp =0.1 the insolvability can become quite substantial 
(0-7 ;$ J? ^ 0.9) for some regions in the (W, cos c .m.) space. 
In other regions of the phase space the insolvability vanishes. 
For ffe X p =0.1, TJimaginary has the largest contribution to rj. 



In Fig. 5, the effect of improving the experimental resolu- 
tion by an order of magnitude is clearly visible. Indeed, both 
the ^imaginary an d ^incorrect surfaces are significantly reduced. 
The overall insolvability, however, still occupies a fair frac- 
tion of the phase space. It is also striking that the regions for 
I cos C m . | w 1 and/or W approaching the threshold are highly 
insolvable, both for ff exp =0.1 and ff exp = 0.01. 

Figure 6 presents a sample extraction analysis of the 
transversity amplitudes for both ff eX p = 0.1 and ff exp = 0.01. 
At each of the 125 x 81 = 10 125 kinematical points a single 
simulated dataset was generated and subsequently analyzed. 
The samples represent how the result of a measurement of 
the (W,cos0 c m ) dependence of the asymmetries might look 
like. For ff exp = 0.1, a large fraction of the kinematical phase 
space suffers from imaginary or incorrect solutions. Increas- 
ing the resolution to ff exp = 0.01 substantially improves the 
situation. For ff exp = 0.1, the fraction of incorrect solutions is 
rather small and for ff exp = 0.01 their contribution is almost 
negligible. 

Though having the smaller contribution to the insolvabil- 
ity, incorrect solutions are spurious and should therefore be 
eliminated, if possible. Incorrect solutions arise from simply 
marking the solution with the highest CI value as the correct 
solution. Such a procedure, however, does not stand on sta- 
tistically solid grounds. A more conservative approach would 
consist of imposing a tolerance level on the CI values. This 
would mean that the solution with the highest CI value is only 
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FIG. 5: (Color online) The {C x , O x ,E,F} insolvabilities rj = I7i ma ginai-y + ^incorrect an d TJincorrect as a function of W and cos c . m . for two values 
of the input experimental resolution (o" exp = 0.1 and o" exp = 0.01). 
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FIG. 6: Two random samples of the individual events that give rise to the left plots shown in Fig. 5. The gray, black, and white dots correspond 
with imaginary moduli and/or phases, incorrect phases, and correct phases. 



accepted as the correct solution if its confidence interval is 
greater than or equals the tolerance level. Imposing a tol- 
erance level would not come without a cost, however. The 
downside of this 'filtering' procedure would be a substantial 
increase in the overall insolvability. Indeed, all of the correct 
solutions that do not pass the tolerance level would also be 
rejected and therefore lead to an increase in the insolvability. 



Nevertheless, the question remains whether such a tolerance 
requirement can really suppress the fraction of incorrect solu- 
tions entirely? 

In Table VII, the mean confidence intervals are listed for 
both 1000 correct and 1000 incorrect solutions at random 
kinematical points. While the "correct" CI values are higher 
than the "incorrect" ones, for each (7 e xp value listed, the dif- 
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TABLE VII: The mean confidence intervals CI for 1000 correct and 
1000 incorrect solutions at random kinematical points (W < 2500 
MeV and cos c . m . e [- 1 , 1]). 





a( 


%) 


Correct 


Incorrect 


0.1 

0.01 

0.001 


93.0^-4 

91 5+ 53 
91 7+5.0 


91 R+ 54 

88.1+^0^2 
84.9+ 2 °;0 



ference between both is not at all substantial. Moreover, it 
is seen that a fair fraction of the incorrect solutions have a 
confidence interval of at least 95%. As there is a significant 
overlap between the "correct" and "incorrect" CI values, the 
entire elimination of the incorrect solutions would require a 
confidence level so high that nearly all of the correct solutions 
would get rejected as well, thereby resulting in a practically 
100% insolvability. Even lowering the confidence level so as 
to remove the bulk of the incorrect solutions, would equally 
wipe out the majority of the correct solutions and hence sub- 
stantially increase the insolvability. 

So, the spurious incorrect solutions can never really be 
identified, at least not at achievable experimental resolutions. 
It was shown in the above analysis that even when the confi- 
dence interval of a solution is significantly high, it is simply 
not possible to state whether the solution in question is a cor- 
rect or an incorrect one. This is in stark contrast with the ob- 
servation that for a considerable experimental resolution the 
fraction of incorrect solutions is nearly negligible. 

V. CONCLUSIONS 

In this paper, the issue of extracting complete informa- 
tion about reaction amplitudes from pseudoscalar-meson pho- 
toproduction data is addressed. The merits of employing 
the transversity basis (for finite experimental resolution) have 
been highlighted. Indeed, linear equations connect the moduli 
of the amplitudes to single-polarization observables. Nonlin- 
ear equations connect the relative phases of the amplitudes to 
double-polarization observables, which are less readily avail- 
able. 

The observables of pseudoscalar-meson photoproduction 
have been expanded in the transversity basis. An inconsis- 
tency with the existing literature is discovered. An indepen- 
dent test of the derived expressions is presented. By rotat- 
ing the obtained transversity expressions to both the helicity 
and the CGLN basis, expressions from literature are retrieved. 
Therefore, convincing evidence has been provided that the de- 
rived transversity expansion for the observables is correct. 

The extraction of the moduli r,- of the normalized transver- 
sity amplitudes is a rather straightforward procedure. The au- 
thors have performed this analysis for the yp — > K + A reaction 
with sets of {L,T,P} data from the GRAAL collaboration. 
These data cover the range 1.65 < W < 1.91 GeV. For the 
moduli, imaginary solutions can be obtained at some isolated 



kinematics due to finite experimental resolution. In the per- 
formed analysis of the GRAAL data, less than 6.5% of the 
extracted r, are complex. Upon improving the experimental 
resolution, the fraction of complex % can be reduced. It is 
found that the (W,cos c .m.) dependence of the extracted?, can 
be nicely reproduced by the RPR-201 1 model. 

The formalism of Wen-Tai Chiang and Tabakin for solving 
complete sets has been extended so as to provide consistent 
estimates for the independent relative phases for data with fi- 
nite error bars. As a check of this formalism, it was applied 
to Monte Carlo simulations of the complete set {C x , O x ,E,F} 
for the reaction yp — » K + A. The simulations are generated 
with the RPR-201 1 model. 

Estimating the independent phases is far more challenging 
than determining the moduli. In order to quantify the phase- 
related issues, the insolvability of the transversity amplitudes 
has been introduced. The insolvability is a measure for the 
fraction of complete measurements which does not result in a 
successful determination of the transversity amplitudes. The 
insolvability receives contributions from both imaginary and 
incorrect solutions. It was observed that the fraction of imagi- 
nary solutions is much larger for the phases than for the mod- 
uli. However, the amount of imaginary solutions can be re- 
duced by increasing the experimental resolution. The 'incor- 
rect' component of the insolvability is much more trouble- 
some. It originates from a discrete phase ambiguity that can- 
not be resolved for finite uncertainties of the asymmetries. Al- 
though this component is not the dominant one and decreases 
for increasing experimental resolution, it was found that at 
achievable experimental conditions it is impossible to discrim- 
inate between correct and incorrect solutions for the phases 
with statistical significance. 

It remains to be investigated whether the measurement of 
an additional double asymmetry, or multiple ones, could help 
in resolving the phase ambiguities for realistic experimental 
resolutions. 
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Appendix: Other representations of the asymmetries 
1. Helicity representation 

The helicity amplitudes H\ : H2 1 H^. and H\ are defined as 

ffi=*(+|/+|-)r, 
H 2 =r(+\J+\+)t, 
H 3 = r (-\J+\-) t , 

H 4 = R (-\J + \+) T . (A.l) 

Here, |±)7- and |±)^ represent the target and recoil helicity 
eigenstates, respectively, and J + is defined in Eq. (9). As the 
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TABLE VIII: The expressions for the single and double asymmetries 
in the normalized helicity basis. 



Helicity representation 



£ 


-Rz(h\h* A -h2h\) 




T 


—\m(h\h2 + hT,h\) 




P 


— Im(/ii/!j + h 2 h^) 




c* 


-Re(h l h^+h 2 hl) 




Q 


4(N 2 + N 2 -N 2 - 


N 2 ) 


o x > 


+Im(/ij/!2 — hsh^) 




Z ' 


— Im(/ii/zJ — h 2 h^) 




E 


-l ( N 2 - N 2 + N 2 - 


N 2 ) 


F 


— Re(/ii^2 + 




G 


+lm(h l hl + h 2 h^) 




H 


+Im(hih^ — h 2 h^) 




T* 
7> 


+Re(h l h* 4 + h 2 h^) 
+Re(/i]/!2 — ^3^14) 




L v , 


-Re{h 1 h^-h 2 hl) 




L z i 


-1 (N 2 -N 2 -N 2 ^ 


-N 2 ) 



The 'normalized' helicity amplitudes hi 

V2H, 



hi 



V\Hi\ 2 + \H 2 \ 2 + \H^ + \H 4 \^ 



(A.7) 



can be written in terms of the normalized transversity ampli- 
tudes a, of Eq. (17) 



(A.8) 



Hereby, use is made of the unitarity of U. From the above 
relations and the expressions of Table I, one obtains the helic- 
ity representation (Table VIII) of the polarization observables. 
The {Ay, A,/} (A e {C,0,T,L}) are related to the {A. T ,AJ 
through 



/ cos 9 C .111. 

— sin 9 C 

in. 

I sin c m . cos c . m . 



(A.9) 



The expressions of Table VIII coincide with those in Ref. [14]. 



2. CGLN representation 



target (recoil) momentum pj (pf>) is directed along the nega- 
tive z-axis (z'-axis), it follows that 



|±>r = |±}- s = «±(pr,M), 

\±)r = \±)-z! = u±(pR,n~ 9 c . m .,n) 



(A.2) 



From Eqs. (2), (3), (4), and (A.2) one deduces the following 
relations 



i±>*=± 



Th 1 ■ 

t.m./2| 



V2 



± 



,±i6U/2 



-F) 



(A.3) 



and in combination with the definition of /_ in Eq. (9), one 
can readily show that the following properties hold 

H 1 = r {+\J + \-) t = +r{-\J-\+)t, 
H 2 = r(+\J+\+)t = -r(-\J-\-)t, 
H 3 =r{-\J+\-)t = -r{+\J-\+)t, 
H 4 =r(-\J+\+)t = +r(+\J-\-)t- (A.4) 

From Eqs. (10), (9), (A.l), and (A.3), one then obtains 

1 



with U a unitary matrix (U^U = UU^ = l) 



(A.5) 



U = 



,i8c.m./2 

~2 



/ 1 


e -iOc.m. 


-1 - e - i9 »M.\ 


i 




i -ie~ i6cm - 


i 




-i ie~ l9cm - 


V-i 


g—i8c.m. 


-1 - e - w ™ J 



(A.6) 



The Chew-Goldberger-Low-Nambu (CGLN) amplitudes Fj 
are, for example, defined in Eq. (8) of Ref. [12]. In what fol- 
lows a connection is established between the the CGLN and 
the transversity amplitudes, as this might be useful for future 
analyses. From the definition of the Dirac spinor (2), the re- 
duced CGLN amplitudes = ^/poFi can be related to the a, 
of Eq. (17). Here, po is the density-of-states factor, defined in 
Eq. (5) of Ref. [12]. The following relation holds 



fi = Vij<*j, 
with V a non-unitary matrix, given by 



V 



i sin C m i sin 9, 



sin 2 9 C 



^c.m. 



cm. 
iflkm. 



V 



1 



-1 






giBcm. 
-1 



(A. 10) 



\ 



- 1 J 

(A.ll) 



This means that the CGLN basis is not orthogonal. By substi- 
tuting a,- = (V~ l )ijfj in Table I, with 










\ 


i 


1 -g'^m 








2 


-1 e- i6 <™ 


— j'sin0 cm . 


_j'g-i8c.m. sin0 c . m . 




\_l e ie c . m . 


/sin0 cm . 


igiSc.m. sm c m y 



v- x = 



(A.12) 

one retrieves the CGLN expansion of the asymmetries as is 
listed in Eqs. (58b-p) of Ref. [12]. 

As V °c sin~ 2 C m ., it is divergent at cos c .m. =±1. In these 
two cases, the CGLN amplitudes cannot be expressed in terms 
of the transversity ones. Conversely, for cos0 c m . = ±1 the 
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a,- can be expanded in the /, basis, though. By substituting 
cos 9 c .m. = ±1 in a,- = {V~ l )ijfj, it is found that 

a\ = -a 2 =a 3 =a 4 = -r(/l T/2)- (A.13) 

The normalization condition (18) then leads to |a t -| = j for 
cos C m . = ±1. As a consequence, by invoking Table I, all the 



single and double asymmetries can be quantified at these two 
extreme angles. At cos0 c m . = ±1 one has that C z = E = 1 
and L z = — I, while all other asymmetries vanish. This result 
holds for any possible value of W. It is quite remarkable that 
this interesting general result can be derived by altering the 
representation of the amplitudes. 
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